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We discuss the mechanism through which classicahzation may occur during the collapse of a 
spherical field configuration modelled as a wavepacket. We demonstrate that the phenomenon 
is associated with the dynamical change of the equation of motion from a second-order partial 
differential equation of hyperbolic to one of elliptic type. Within this approach, we rederive the 
known expression for the classicalization radius. We also find indications that classicalization is 
associated with the absence of wave propagation at distances below the classicalization radius and 
the generation of shock fronts. The full quantitative picture can be obtained only through the 
numerical integration of a partial differential equation of mixed type. 
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The classicalization scenario advocates that the high-energy behavior of certain classes of seemingly simple 
scalar field theories can be nontirivial. Typically, such theories are described by Lagrangians which include non- 
renormalizable higher-derivative terms. It is possible that scattering in the context of such theories can take place at 
^^"^ , distances much larger than the typical length scale associated with the couplings multiplying the non-renormalizable 

■ terms. The scattering scale , which is termed classicalization radius and determines the cross-section, depends on 

the center of mass energy in a nontrivial way. The conjecture put forward in ref. is that distances much smaller 

^ • than r* cannot be probed, as all significant scattering takes place already at r*. Thus, the loss of unitarity expected 
Oh when probing length scales ~ L*, as a result of the presence of the non-renormalizable terms, never occurs. 
^ ■ The classicalization picture relies heavily on the analogy with black-hole formation at ultra-Planckian energies. 
The high-energy scale is the analogue of the Planck scale, while the classicalization radius the analogue of the 
Schwarzschild radius. The classicalon is a configuration of the scalar field that mimicks the behavior of the black hole. 
\ Both objects are one-parameter solutions of the equations of motion with a singularity at the origin. In the case of the 
^ ■ black hole the singularity is hidden behind the horizon. For the classicalon, it is usually attributed to a source term 
[ generated by the self-sourcing of the field in situations in which its energy is concentrated in small regions of space. 
T— I ■ In specific examples, the field configuration is split into a part 4>o{t, x), that corresponds to a solution of the equation 
' of motion in the absence of the non-renormalizable terms, and a perturbation <pi{t, x). Perturbation theory generates 

• an equation of motion for (f>i , which includes a source term resulting from the variation of the non-renormalizable 
' terms evaluated for ^q. 

• It seems clear that the perturbative arguments need to be reinforced by a non-perturbative analysis. As the essence 
" ^ I [ of the classicalization picture is that only classical physics is relevant for the theories in question, the necessary 

J> ■ procedure is obvious: One needs to solve the classical equation of motion of the full theory for the scattering problem. 
', The purpose of this letter is to study the form of this equation for the case of a collapsing spherical wavepacket 
■ discussed in ref. 0, 0] ■ Even though an accurate numerical solution is a difficult technical task that we postpone for 
d \ the future, the general properties of the equation can be established more easily. In particular, the emergence of the 
classicalization radius r* , the formation of shock fronts and the occurence of scattering already at this scale, as well 
as the suppression of wave propagation at length scales smaller than r, , seem plausible predictions of the equation of 
motion. Our analysis does not rely on the emergence of a configuration identified with the classicalon. Nevertheless 
it supports the picture of scattering at length scales much larger than the short-distance scale L*. 



EQUATION OF MOTION 

We consider the theory of a scalar field x) with derivative self-interactions. The Lagrangian density is 



(1) 



with the Minkowski metric given by 77^1, = diag(l, —1, —1, —1). We allow for both signs of the higher-derivative term 
by considering the values 6 = ±1. The equation of motion of the field (f) is 



1 - 5Ll {d. 



= 0. 



(2) 
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We are interested in a configuration described by a collapsing spherical wavepacket represented by a Gaussian of 
width a initially centered around a radius tq. It has the form 



/ , A 
4'o{t,r) = — exp 
r 



[r + t-rof 



(3) 



This configuration is an exact solution of the equation of motion in the absence of the higher-derivative term (L» = 0). 
It is also an approximate solution for nonzero L* and sufficiently large values of tq. 
When expressed in spherical coordinates, eq. 1^ assumes the form 

(1 - 3A(/.2 + X(bl) ^tt - (1 - + 3X(bl) <prr + 4A0,0t q^tr = ^ (l - A^^ + A02) , (4) 

where A = SL'^ and subscripts denote partial derivatives. The above equation can be expressed as a conservation law, 
in the form 

dt [cf>t (1 - A02 ^ ^^2)] _ }_Q^ ^^2^^ _ ^^2 ^ ^^2^] ^ Q (5) 

The initial conditions for its solution are of the Cauchy type: 

0(O,r) = 0o(O,r), dtmr)=dtMO,r), (6) 

where we assume that is much larger than any other physical scale. 

In refs. 0, [1] eq. @ is solved by writing (f>{t,r) = (j)Q{t,r) + (j)i(t,r), with 0o given by eq. ([3]). The correction 
is treated as a perturbation that vanishes for — 0. It is shown that (/>i becomes comparable to (jjo Ett the 
classicalization radius 



(7) 

The subsequent evolution of the field cannot be derived reliably within perturbation theory. Moreover, it is not clear 
whether the initial configuration, despite its deformation, continues its propagation towards smaller values of r, or 
whether there is significant scattering towards large r. 

The properties of eq. (U) become more transparent if it is written in the form of a quasi-linear second-order partial 
differential equation: 

A{4>t, 4>r) 4>tt + B{4>t, (t)r) 4>tr + C{(l)t,4>r) (t>rr = 'D{4)t,(j)r, r) , (8) 

with 

= l-3A02-fA02 (9) 

B{(f,t,^r) = 4XMr (10) 

C{^t,A) = - (1 - A02 + 3A02) (11) 
V{cf>„cf>r,r) = ^ (1 - A02 + A,^2) ^ (12) 

The type of this partial differential equation is determined by the discriminant 

A = B^ - 4:AC = 12 



Q - -I- A02 j (1 - + A02) . (13) 



For A > the equation is hyperbolic, for A = parabolic, while for A < elliptic. Under an appropriate coordinate 
transformation {t, r) — >■ (a, (3), a hyperbolic partial differential equation can be written in the form ipaa — 4>i3[3 + ... = 0, 
while an elliptic one as (j)aa +4'pi3 + •■• — Oj with the dots denoting lower-derivative terms. Hyperbolic equations have 
wave-like solutions, while elliptic ones do not support propagating solutions. It is clear that, depending on the form 
of (j){t, r), eq. dH) can be of any type. Moreover, it is possible that it may change type dynamically during the time 
evolution of an initial configuration. Equations that have this property are characterized as of mixed type and are 
notoriously difficult to integrate. 
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FIG. 1: The characteristics of the field equation ([8|. The functions A{4>t,4'r), S(0t,<^r), C{(l)t,4>r) are evaluated for the 
configuration 4>o{r,i) of eq. ((3]) , with a = f , A = 10, ro = 15, in units such that L» = 1, A = 1. The equation is of the elliptic 
type within the shaded (red) area. The panel (b) is a magnification of the area denoted by a square in panel (a). 



For large r we have A{4>ot, (por) — C{<pot, 0Or) — 1 and B{(/)ot, 4>or) — 0, where A, B, C have been evaluated for the 
configuration (/)o(t, r) of eq. ([3]). In this limit, A ~ 4 and eq. ([8]) becomes the wave equation in spherical coordinates. 
The configuration ([3]) is an approximate solution, as it is a superposition of spherical waves. It is possible however, 
that the subsequent evolution of (j>{t, r) may result in the change of the type of eq. ([S]). We demonstrate this possiblity 
by depicting in fig. [ijthe region in which A < 0, when A, B, C are evaluated for the configuration (jSj with a = 1, 
A — 10, To — 15. We consider the case with 6 = 1 and work in units such that = 1 and A = 1. The region in which 
eq. ([5]) becomes elliptic corresponds to the shaded (red) area on the {t, r) plane. It must be emphasized that eq. @ 
is not a solution of eq. (|H]) for late times and small values of r, so that the change of equation type must be viewed 
only as a possibility. The derivation of an accurate (numerical) solution of eq. (|S]) with initial conditions given by ^ 
is necessary in order to obtain a definite confirmation of this conclusion. 

It is interesting to estimate the distance r at which eq. (|8]) changes type if A, B, C are evaluated for the configuration 
([3|). The discriminant vanishes for 

= -k'XA'^e-'^'^ (1 + 4fe) , (14) 

where r = r/a, fo = fo/o.^ ^ = ^/ci'^, e = f + t — fg and k takes the values 1 or 3 (corresponding to the inner and 
outer curves bounding the shaded (red) region in fig. [l}. The relevant values of the parameter e are of order 1, as 
the Gaussian falls off very rapidly outside this range. We are interested in the region f 3> 1, because we expect the 
classicalization radius to be much larger than the width of the Gaussian. A solution can be obtained for e ~ — 1 and 
f|e| >• 1. We obtain 

h [XA^j , (15) 

an expression which coincides with the definition of the classicalization radius r* in eq. ([7]). It is clear that can be 
interpreted as the radius at which the partial differential equation ([8]) is expected to switch from hyperbolic to elliptic. 
The change does not occur simultaneously at all values of r. For example, the center of the Gaussian corresponds to 
e ~ 0, so that eq. ([T4l) is never satisfied for A > 0. This indicates that the equation changes type first at the front of 
the wavepacket, i.e. for values of r smaller than the one corresponding to the center of the Gaussian. 

The analysis can be repeated for 5 = — 1, which implies A < and corresponds to a higher-derivative term with 
positive sign in the Lagrangian ([1]). The change of type of the equation of motion also occurs in this case. In fig. [2]we 
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FIG. 2: The characteristics of the field equation ^ for a = 1, A = 10, ro = 15 and A — —1. The panel (b) is a magnification 
of the area denoted by a square in panel (a) . 

depict the region in which A < 0, when A, B, C are evaluated for the configuration ^ with a = 1, A = 10, vq = 15 
and A = — 1. The region in which eq. (jS]) becomes elliptic corresponds to the shaded (red) area on the {t,r) plane. 
This area is similar to the one depicted in fig. [l] but slightly displaced. The analysis of the previous paragraph, 
leading to the estimate of the classicalization radius, can be repeated for A < 0. The result is 

f,^(|A|A2)'^', (16) 

which again reproduces eq. ([7]). The main difference is that the type change now occurs first at e ~ 1, i.e. at a point 
at the rear of the wavepacket. Also, in this case there is a solution of eq. (IT4l) for e = 0, given by fc ~ (|A| A^)-'^/'*. For 
f» 3> 1 we have f» > f^. 

CHARACTERISTICS 

In order to obtain more intuition on the nature of the solutions of eq. (|8]), we examine its characteristics when 
A{(f>t,(f>r), B{(j)t,(f>r), C{(j)n4'r) arc evaluated for the configuration (f)i^{r,t) of eq. ([3|). The characteristics are solutions 
of the equation 

Clearly, they can be defined only for A > 0, and they do not exist in the region where the equation is elliptic. The 
two signs correspond to two families of characteristics. For large r and early times we have y^~C~ 1, S~0. Eq. 

becomes the wave equation, which has the configuration (j)o{t,r) as an approximate solution. The characteristics 
correspond to the lines r -\- 1 = ci, r — t = C2. The configuration (poit, r) is projected only on the first family, as it is 
composed from ingoing spherical waves, which are functions of r + 

The deformation of the characteristics in the vicinity of the region in which the equation changes type is depicted 
in fig. [T]for a model with a = 1, A = 10, ro = 15. We consider the case 5 = 1 and use units such that L* = 1, A = 1. 
The two families correspond to solutions of eq. (fTTl) expressed in the form u{t, r) = ci and v{t, r) = C2. Each value of 
ci or C2 determines a curve on the (i, r) plane, along which initial disturbances propagate. We choose the first family 
(which we term ingoing) as the one that reduces to u(t, r) = r + 1 at initial times and large r, while the second one 
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(termed outgoing) reduces to v{t, r) = r — t. The dash-dotted (black) hnes in fig. [T] depict ingoing, while the solid 
(green) lines outgoing characteristics. 

The initial configuration (f>Q(t,r) is projected only on the first family and within a limited range of r. For t ~ 
the initial Gaussian is centered at tq = 15 and has a width a — 1. The initial disturbance is localized roughly in the 
region 14 ^ r ^ 16. At later times it propagates along the ingoing characteristics that start within this range. We 
have indicated its position by emphasizing the ingoing characteristics starting at r = 14 and r = 16, which roughly 
bound the disturbance. The ingoing characteristics remain linear to a good approximation, even in the vicinity of the 
shaded (red) region. They terminate on the boundary of this surface. The characteristics starting at r < 15, which 
correspond to points at the front of the wavepacket, reach the shaded region before the center of the Gaussian, which 
starts at r — 15. 

The outgoing characteristics are deformed strongly in the vicinity of the shaded area. Even though they are 
equidistant near the r-axis, they come very close, almost overlapping, near the shaded area. This indicates that 
distinct points of an outgoing configuration tend to merge when they pass through this area. Another striking feature 
is that the "velocity" dr/dt diverges at certain points along the characteristics. All these points are located on the 
lines determined by the condition ^ = 0, which are depicted by dashed thick lines in fig. [TJ The form of the outgoing 
characteristics is visible more clearly in the enlarged region (part (b) of fig. [1} which corresponds to the square in 
part (a). This form is consistent with the appearance of outgoing shock fronts. We have depicted only characteristics 
emanating from the left of the shaded (red) region in order to demonstrate their convergence. There are also outgoing 
characteristics starting on the boundary of the (red) region and covering the right part of the diagram. We do not 
depict those for clarity of the picture. 

The form of the characteristics for a model with a = 1, A = 10, ro = 15 and A = — 1 is depicted in fig. [2j The 
ingoing characteristics remain linear to a good approximation, similarly to the case with A = 1. They terminate on the 
boundary of the shaded (red) surface. The characteristics starting at r > 15, which correspond to points at the rear 
of the wavepacket, reach the shaded region before the center of the Gaussian, which starts at r = 15. The outgoing 
characteristics are again deformed strongly in the vicinity of the shaded area. Even though they are equidistant near 
the r-axis, they approach each other near the shaded area. Distinct points of an outgoing configuration tend to merge 
when they pass through this area. Again this suggests the formation of a shock front. However, the various points 
subsequently move away from each other and the characteristics become dilute at large t. Another striking feature is 
that the "velocity" dr/dt vanishes at the points along the characteristics where C — 0. These are depicted by dashed 
thick lines in fig. [2] 

CONCLUSIONS 

Our study of the partial differential equation ([8|) has been based on evaluating the coefficient functions A{(l)t ■, 4'r) , 
B{4it,4>r), C{<j)t,<j)r) for the spherical wavepacket (j)(){r^t) given by eq. ([3]). This configuration is an approximate 
solution of the exact eq. for early times and large r. However, it is not expected to remain so when the center 
of the wavepacket approaches the region near r = 0. The main prediction of the classicalization scenario is that the 
deviations become large already at the classicalization radius given by eq. ([7]), and not at the fundamental scale 
i*. We have analyzed this possibility by examining the nature of the solutions of eq. (|8|). 

Our main observation is that eq. ([8]) is a partial differential equation of mixed type. For the initial conditions 
^ the equation is hyperbolic in the part of the (t, r) plane where the field is nonvanishing initially. At later times 
the field fiuctuation moves into a region where the equation becomes elliptic. Wave propagation is not supported by 
equations of elliptic type, and is probably replaced by exponential decay. On the other hand, it seems also plausible 
that the inability of the fluctuation to propagate towards smaller values of r will be accompagnied by, at least partial, 
reflection. The situation is complicated by the indications for outgoing shock-front formation that we discussed in the 
previous section. We expect that these outgoing shock fronts will form flrst at the front or rear of the wavepacket, 
depending on the sign of the higher derivative term in the Lagrangian ([1]). Subsequently, they will interact with the 
remaining part of the ingoing wavepacket. Clearly, our current analysis cannot provide any clues on the details of 
this process. We also mention at this point that the numerical treatment of the formation of shock fronts is a very 
difficult problem, for which a general framework has not been developed yet. Instead, the analysis is carried out on 
a case by case basis. 

Despite neglecting the effect of backreaction in the results of the previous sections, we may attempt to draw some 
quantitative conclusions. The ingoing characteristics in flg. [T] remain largely unaffected up to the point where they 
terminate on the surface A = 0. As a result, we expect that the reflection of (part of) the wavepacket will take place 
within a short range of values of r. The relevant scale is given by the classicalization radius of eq. ([7]). Also, the 



divergence of the speed of transmission observed at certain points on the outgoing characteristics is an artefact of 
our approximate treatment of eq. ([8]). In an exact solution, causality would prevent this divergence. However, the 
merging of characteristics, typical of the appearance of shock fronts, is expected to persist. 

As a final comment, we mention that the approach that we outlined is applicable also to cases in which the 
phenomenon of classicalization is not expected to occur. For example, for planar configurations moving along the 
z-axis the equation of motion becomes 



(1 - 3A</)? + Ml) - (1 - X<Pt + 3A0^) (j)^, + iX(j),(j)t <pt. = 0. 



(18) 



An initial field configuration of the form 



io(i, z) = Aexp 



Aexp 



{z-t + Zo) 



(19) 



describes two localized planar wavepackets moving in opposite directions. It is an approximate solution of eq. (jlSp 
as long as there is no significant overlap of the two Gaussians. The discriminant of the partial differential equation 
(UHl) is 
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Q - Ml + Ml^ (1 - Ml + Ml) . 



(20) 



Evaluating it for the configuration ([19]) gives A ~ 4 (the equation is hyperbolic), apart from the region of significant 
overlap of the wavepackets. This is is a region of width ~ a around the origin, appearing at times t zq. The 
classicalization phenomenon is not expected to occur in this case, in agreement with the analysis of ref. 
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